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ABSTRACT 

We present a study of the individual super star clusters (SSCs) in the low-metallicity galaxy SBS 
0335-052 using new near-infrared and archival optical Hubble Space Telescope observations. The 
physical properties of the SSCs are derived from fitting model spectral energy distributions (SEDs) 
to the optical photometry, as well as from the Ha and Paa nebular emission. Among the clusters, 
we find a significant age spread that is correlated with position in the galaxy, suggesting successive 
cluster formation occurred in SBS 0335-052 triggered by a large-scale disturbance traveling through the 
galaxy at a speed of ~ 35 km s _1 . The SSCs exhibit /-band (~ 0.8 /xm) and near-IR (~ 1.6 — 2.1 jim) 
excesses with respect to model SEDs fit to the optical data. We hypothesize that the /-band excess is 
dominated by a photoluminescent process known as Extended Red Emission; however, this mechanism 
cannot account for the excesses observed at longer near-IR wavelengths. From the cluster SEDs and 
colors, we find that the primary origin of the near-IR excess observed in the youngest SSCs (< 3 Myr) 
is hot dust emission, while evolved red supergiants dominate the near-IR light in the older (> 7 Myr) 
clusters. We also find evidence for a porous and clumpy interstellar medium (ISM) surrounding the 
youngest, embedded SSCs: the ionized gas emission underpredicts the expected ionizing luminosities 
from the optical stellar continuum, suggesting ionizing photons are leaking out of the immediate 
vicinity of the clusters before ionizing hydrogen. The corrected, intrinsic ionizing luminosities of the 
two SSCs younger than ~ 3 Myr are each ~ 5 x 10 52 s _1 , which is equivalent to each cluster hosting 
- 5000 07.5 V stars. The inferred masses of these SSCs are - 1O 6 M . 

Subject headings: galaxies: dwarf - galaxies: individual (SBS 0335-052) - galaxies: starburst - galaxies: 
star clusters 



1. INTRODUCTION 

Ancient globular clusters are among the oldest objects 
known, almost as old as the universe itself with ages es- 
timated in excess of 12 Gyr (Vandenberg et al. 1996; 
Freeman & Bland-Hawthorn 2002). As such, these clus- 
ters are valuable relics of the earlier universe, when vi- 
olent and intense star formation was common. Indeed, 
globular clusters are ubiquitous around massive galaxies 
today (Harris 1991; Brodie & Strader 2006), and they 
must have been formed prodigiously in the primordial 
universe given that > 90% may have been subsequently 
destroyed (Fall & Zhang 2001; Whitmore et al. 2007). 

However, the conditions required for the creation of 
globular clusters have puzzled astronomers for decades. 
For many years the prevailing belief was that globular 
clusters were simply formed by the gravitational collapse 
of density inhomogeneities in the early universe (Peebles 
& Dicke 1968; Fall & Rees 1985), and little was known 
about their early evolution. However, in the mid 1990's, 
observations using the Hubble Space Telescope discovered 
extremely young, massive, and compact clusters in the 
local universe, the so-called "super star clusters" (SSCs) 
that have properties consistent with those expected for 
adolescent globular clusters (e.g. Whitmore 2003). The 
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discovery of SSCs precipitated a major shift in our under- 
standing of the conditions required for globular cluster 
formation. 

Although the evidence connecting young SSCs and an- 
cient globular clusters is compelling, a critical issue re- 
mains: modern day SSCs and ancient globular clusters 
were formed in environments with very different metallic- 
ities. We have not yet observationally constrained what 
effect low metal abundances had on the formation of 
massive star clusters. Much effort has gone into under- 
standing the formation of the first stars in the universe 
(e.g. those stars formed out of truly primordial material, 
Bromm & Loeb 2003; Santoro & Shull 2006; Tumlinson 
2007), and it is clear that metallicity has a critical role 
in this regime. Even moderately low metal abundances 
may affect star formation in a number of ways, including 
dust formation, cooling and pressure, and the hardness 
of the resulting stellar spectra. 

In an effort to investigate massive star cluster for- 
mation in an environment similar to that which might 
be found in primordial galaxies during the time ancient 
globular clusters were prolifically formed, we have tar- 
geted the galaxy SBS 0335-052 (E) for a detailed multi- 
wavelength study. SBS 0335-052 is a remarkable blue 
compact dwarf galaxy at a distance of 54 Mpc (NED 1 ) 
that is well-known for its extremely low oxygen abun- 
dance of ~ 1/40Z© (Izotov et al. 1990; Melnick et al. 

1 NASA/IPAC Extragalactic Database 
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Fig. 1. — HST images of the observations used in this work. 



NIC2 Pact 



ACS/WFC FR656N 



* 

* * # 



f 



ACS/HRC F330W 



NIC2 F1 60W 
ACS/WFC F550M 

m 

* 

ACS/HRC F220W 



1992; Izotov et al. 2001). Unlike I Zw 18 and SBS 0335- 
052W, which have slightly lower metallicities, SBS 0335- 
052 is undergoing a vigorous starburst with a star for- 
mation rate > 1 M yr _1 and hosts extremely massive 
young clusters, making it ideal for this study (Thuan et 
al. 1997, 1999; Hunt et al. 2001; Dale et al. 2001; Plante 
& Sauvage 2002; Hunt et al. 2004; Houck et al. 2004). 

Here we present new near-infrared and archival optical 
high-resolution observations of SBS 0335-052 with the 
goal of studying SSC formation and evolution at low 
metallicity. Near-IR observations are particularly impor- 
tant for probing the dusty birth environments of embed- 
ded massive star clusters, and the wavelength range of 
~ 1 — 2 /im is well-suited to assessing a cluster's evolu- 
tionary state because it samples the spectral energy dis- 
tribution (SED) at the nexus of emission between dust 
and stellar light. 

The HST observations and photometry are presented 
in §2 and §3. In §4, we discuss various properties of the 
SSCs. Specifically, we derive the physical properties of 
the SSCs from fitting model SEDs to the optical pho- 
tometry, as well as from the Ha, Paa, and thermal radio 
nebular emission (Johnson, Hunt, & Reines in prep). In 



addition, we present evidence for a porous and clumpy 
ISM in the youngest embedded clusters in SBS 0335-052 
which can account for the apparently discrepant extinc- 
tion estimates found in the literature. We also provide 
explanations for the observed /-band (~ 0.8 /im) and 
near-IR (~ 1.6 — 2.1 /im) excesses. Evidence for suc- 
cessive cluster formation is given in §5, and finally, our 
conclusions are summarized in §6. 

2. HST OBSERVATIONS 

Mult i- wavelength imaging of SBS 0335-052 has been 
obtained with the Hubble Space Telescope 2 . New high- 
resolution near-IR observations are presented along with 
archival UV and optical images of the galaxy. A sum- 
mary of the observations used in this work is given in 
Table 1 and the images are shown in Figure 1. The 
SSCs (1-6) are identified in Figure 2 using the notation 
of Thuan et al. (1997). 

2 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with program #10894, 10575, & 5408. 
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TABLE 1 
HST Observations of SBS 0335-052 



Filter 


Instrument 


Description 


Exp. Time (s) 




New Near-IR Observations a 


F160W 


NIC 2 


1.6 /im, ~H 


2431 


F187N 


NIC 2 


Pao; continuum 


4862 


F190N 


NIC 2 


Pa (redshifted) 


4606 


F205W 


NIC 2 


2.1 /im, ~K 


2175 


Archival UV and Optical Observations 13 


F220W 


ACS/HRC 


Near-UV 


1660 


F330W 


ACS/HRC 


HRC U 


800 


F435W 


ACS/WFC 


Johnson B 


680 


F550M 


ACS/WFC 


Narrow V 


430 


FR656N 


ACS/WFC 


Ra 


680 


F791W 


WFPC2 


I 


4400 



a The NIC 2 observations are associated with Proposal ID 
10894 (this work). 

b The ACS and WFPC2 observations are associated with 
Proposal IDs 10575 (PI G. Ostlin) and 5408 (PI T. Thuan), 
respectively. 



2.1. Near-IR Imaging with NICMOS 

High-resolution near-IR observations of SBS 0335-052 
were obtained with the Near Infrared Camera and Mul- 
tiobject Spectrometer (NICMOS) on 2006 September 28 
and 2006 December 25-26. Using the NIC 2 Camera, 
broad-band images at 1.6 /im (~ H) and 2.1 fim (~ K) 
were obtained through the F160W and F205W filters, 
respectively, and narrow-band Paschen a line and con- 
tinuum images were obtained through the F190N and 
F187N filters, respectively 3 . A spiral dither pattern was 
implemented to allow improvement of the spatial reso- 
lution of the observations by better sampling the point 
spread function (PSF) of the NIC 2 Camera. In addi- 
tion, off-source observations were obtained in the F187N, 
F190N, and F205W filters in order to mitigate the effects 
of the thermal background. 

The raw data were processed with the calnica NIC- 
MOS calibration pipeline. The calibrated NICMOS data 
were further reduced using a combination of custom IDL 
programs and the Drizzle software (Fruchter & Hook 
2002) available in the STSDAS 4 DITHER package in 
IRAF 5 . The residual flat imprint and "pedestal effect" 
that can plague NICMOS observations were effectively 
removed by subtracting off-source "sky" images from on- 
source images for each dither position 6 and then sub- 
tracting the quadrant-dependent bias levels. The result- 
ing images (8 for each of the F160W and F205W filters 
and 16 for each of the F187N and F190N filters) were 
then combined using the Drizzle software to remove cos- 

3 At the redshift of SBS 0335-052, the Paa line falls in the F190N 
filter and the F187N filter is used to obtain the continuum. 

4 The Space Telescope Science Data Analysis System 

5 The Image Reduction and Analysis Facility is distributed by 
the National Optical Astronomical Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy 
(AURA) Inc., under cooperative agreement with the National Sci- 
ence Foundation (NSF). 

6 This step was not applied to the F160W images since thermal 
fluctuations are minimal at 1.6/i,m and we did not chop to sky for 
these observations. 




Fig. 2. NICMOS F160W (~ if-band) image with the SSCs 
identified by Thuan et al. (1997). 

mic rays and improve the spatial resolution of the images. 
The drizzle parameter scale was set to 0.5 to reduce the 
native NIC 2 plate scale of ~0"075 pixel -1 to ~0'/0375 
pixel -1 . The final drizzled images are essentially diffrac- 
tion limited, with spatial resolutions in the range ~(//15 
(1.6 /im) to ~0'/2 (2.1 /im). Figure 3 shows a three- 
color NICMOS image: F160W (blue), F205W (green), 
and Paa (red). 

2.2. Archival UV and Optical Images 

Optical and UV observations of SBS 0335-052 have 
been obtained with the Advanced Camera for Surveys 
(ACS) and we have retrieved the pipeline-produced cal- 
ibrated and drizzled images from the archive. Images 
through the F220W (~ Near-UV) and F330W (~ U) 
filters were obtained with the High Resolution Channel 
(HRC) and images through the F435W (~ £), F550M 
(~ V), and FR656N (Ha) filters were obtained with the 
Wide Field Channel (WFC). The HRC and WFC have 
plate scales of ~0 / /027 pixel -1 and ~0 / /05 pixel -1 , re- 
spectively. The HRC UV images have a resolution of 
~0'/06 and the WFC optical images have a resolution of 
~0 / /125. SSC 3 is actually two sources at the resolution 
of the ACS/HRC. 

Two images of SBS 0335-052 through the F791W (~ I) 
filter were obtained with the Wide Field and Plane- 
tary Camera 2 (WFPC2) and we combined the pipeline- 
produced calibrated images from the archive and rejected 
cosmic rays. The galaxy was placed on the WF3 chip 
which has a plate scale of ~0'/l pixel -1 . The resolution 
of the F791W image is ~0'/2. 

3. PHOTOMETRY OF THE SUPER STAR CLUSTERS IN 
SBS 0335-052 

We performed aperture photometry on the super star 
clusters identified by Thuan et al. (1997) using a custom 
IDL photometry program (SURPHOT) allowing for con- 
sistent apertures and background annuli across multiple 
wavebands. A full description of SURPHOT is given in 
Reines et al. (2008). 
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Fig. 3. — NICMOS three-color image: F160W (~ H) in blue, F205W (~ K) in green, and Pao; (continuum subtracted) in red. 



The SSCs are identified in the F160W image (using 
the notation of Thuan et al. 1997, see Figure 2) and 
the flux density at each wavelength is measured within 
a circular aperture of radius 0'/15. The background 
contribution at each wavelength is determined within 
an annulus of inner and outer radii equal to 0"45 
and 0'/65, respectively. Aperture corrections are de- 
termined from simulated PSFs for each camera-filter 
combination and position on the detector using the 
Tiny Tim software package (v6.3) distributed by the 
Space Telescope Science Institute (Krist & Hook 2004). 
Final aperture-corrected flux densities are obtained 
by multiplying the measured background-subtracted 
flux densities within a 0'/15 radius aperture by 
the following factors: 1.29(F220W), 1.21(F330W), 
1.23(F435W), 1.22(F550M), 1.22(FR656N), 

1.68(F791W), 1.85(F160W), 1.87(F187N), 1.88(F190N), 
and 1.94(F205W). 

Final aperture-corrected flux densities and HST sys- 
tem magnitudes with Vega zeropoints of the SSCs are 
listed in Table 2. Paa fluxes, Ha fluxes, and Ha equiva- 
lent widths of the SSCs are given in Table 3. The width 
of the FR656N (-131 A) filter is such that it would con- 
tain any [Nil] emission in addition to Ha. The [Nil] 



emission in SBS 0335-052, however, is negligible (Izo- 
tov et al. 1990). The Ha equivalent widths are obtained 
by dividing the total (continuum subtracted) Ha flux in 
the FR656N filter by the continuum flux density at the 
redshifted wavelength of Ha (6651 A), which is found 
by interpolating between the F550M and F791W filters. 
We note that the Paa flux measured in the entire line- 
emitting region containing SSCs 1 and 2 is equal to that 
found by Thompson et al. (2006), 1.77 x 10~ 14 erg s" 1 
cm -2 . The Ha flux in this region is 1.66 x 10 -13 erg s _1 
cm -2 . 

4. PROPERTIES OF THE SSCS 

4.1. Spectral Energy Distributions 

We estimate the physical properties of the SSCs in 
SBS 0335-052 by comparing their spectral energy dis- 
tributions (SEDs) to the latest STARBURST99 popula- 
tion synthesis models (v5.1) of Leitherer et al. (1999). 
We have run two simulations using different evolution- 
ary tracks and metallicities: (1) the Geneva tracks with 
high mass loss and Z=0.001 and (2) the Padova tracks 
with AGB stars and Z=0.0004. In both cases, we adopt 
an instantaneous burst of 1O 5 M with a Kroupa IMF 
(Kroupa 2001) and the Pauldrach/Hillier atmospheres. 
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TABLE 2 

Photometry of the SSCs in SBS 0335-052 



Source F220W F330W F435W F550M F791W F160W F187N F205W 







Log 


Flux Density (erg s 


1 cm- 2 I 


L" 1 ) 






SSC 1 


-15.66 


-15.90 


-16.10 


-16.51 


-16.65 


-17.56 


-17.51 


-17.35 


SSC 2 


-15.69 


-15.98 


-16.20 


-16.58 


-16.83 


-17.60 


-17.68 


-17.49 


SSC 3 


-15.86 


-16.28 


-16.53 


-16.86 


-17.11 


-17.88 


-18.09 


-18.08 


SSC 4 


-15.39 


-15.87 


-16.15 


-16.55 


-16.87 


-18.06 


-18.13 


-18.18 


SSC 5 


-15.34 


-15.78 


-15.98 


-16.27 


-16.56 


-17.58 


-17.64 


-17.73 


SSC 6 


-16.31 


-16.71 


-16.91 


-17.18 


-17.41 


-18.37 


-18.54 


-18.46 



HST System Magnitude with Vega Zeropoint 



SSC 1 


18.29 


18.53 


19.78 


20.11 


19.29 


19.16 


18.37 


17.65 


SSC 2 


18.37 


18.73 


20.02 


20.28 


19.75 


19.27 


18.79 


18.01 


SSC 3 


18.79 


19.48 


20.85 


20.98 


20.44 


19.96 


19.82 


19.49 


SSC 4 


17.62 


18.45 


19.90 


20.20 


19.85 


20.41 


19.92 


19.75 


SSC 5 


17.50 


18.24 


19.47 


19.52 


19.07 


19.20 


18.68 


18.61 


SSC 6 


19.91 


20.54 


21.80 


21.79 


21.20 


21.17 


20.95 


20.43 



Note. — We adopt a distance modulus of 33.66 for SBS 0335-052. The uncertainties 
in the flux densities are ~ 10%. 



TABLE 3 
Nebular Emission from the SSCs 



Source 


Paa Flux (xl0~ 15 
erg s _1 cm -2 ) 


Ra Flux (xlO- 15 
erg s _1 cm -2 ) 


Ha Equivalent 
Width (A) 


SSC 1 


7.26(0.78) 


61.17(6.48) 


2300(410) 


SSC 2 


4.65(0.50) 


37.69(4.05) 


1870(330) 


SSC 3 


0.55(0.07) 


6.18(0.78) 


580(110) 


SSC 4 


0.23(0.04) 


1.55(0.56) 


80(30) 


SSC 5 


0.18(0.07) 


1.19(0.96) 


30(20) 


SSC 6 


0.08(0.01) 


0.23(0.13) 


50(30) 



The two simulations use the lowest metallicity evolution- 
ary tracks available in STARBURST99 from the Geneva 
and Padova groups, appropriate for SBS 0335-052. The 
high mass loss tracks are recommended by the Geneva 
group and the inclusion (or exclusion) of AGB stars in 
the Padova tracks is not critical since we are modelling 
very young clusters. We have run the two simulations 
for comparison and as a consistency check. 

Ages, extinctions, and masses of the SSCs are esti- 
mated by comparing a grid of model SEDs (separately 
for each simulation) to the broad-band HST flux densi- 
ties. The grids consist of model SEDs with ages up to 
20 Myr in increments of At = 0.1 Myr and extinctions 
up to Ay = 3 mag in increments of 0.1 mag. The work 
of Reines et al. (2008) on young massive clusters in the 
dwarf starburst galaxy NGC 4449 showed that model 
SEDs modified by a 30 Doradus extinction curve 7 fit 
the clusters' broadband flux densities significantly bet- 
ter than the Galactic extinction curve (Cardelli et al. 
1989) and the starburst obscuration curve (Calzetti et 
al. 2000). Therefore, we first apply the 30 Doradus ex- 
tinction curve with a given Ay to the model SEDs. Next, 
we simulate Galactic foreground extinction towards SBS 

7 The 30 Doradus extinction curve is adopted from Misselt et al. 
(1999, Table 3) and Fitzpatrick (1985, Table 6), using the param- 
eterization given by Fitzpatrick Sz Massa (1990). 



0335-052 (E(B - V) = 0.047 mag, Schlegel et al. 1998) 
using the extinction curve of Cardelli et al. (1989). Each 
model SED, of a given age and Ay, is then convolved 
with the HST total throughput curves of the appropriate 
filters before it is compared to the measured flux densi- 
ties. The best-fit model SED is determined by minimiz- 
ing a goodness-of-fit parameter, a, equal to the standard 
deviation of the logarithmic residuals of the observed and 
model flux densities. A mass estimate is obtained by 
scaling the model mass by the mean logarithmic offset 
between the observed and best-fit model flux densities. 

We include the F220W, F330W, F435W, and F550M 
flux densities in the SED fitting. The strongest Balmer 
lines, Ka and H/3, are not contained in this filter set, al- 
though the F435W filter does include H7 and US. Since 
Reines et al. (2008) revealed excesses in the /- and H- 
band flux densities of the young clusters in NGC 4449 
compared to model SEDs, we did not include data at 
wavelengths longer than the F550M filter (~ V) in the 
fitting process here. The observed and best-fitting model 
SEDs of the SSCs in SBS 0335-052 are shown in Figure 
4. The grids of a are also shown to give a sense of the un- 
certainty in the estimates of ages and extinctions. Black 
indicates lower values of sigma (better fits) and red indi- 
cates higher values of sigma (worse fits). The best-fitting 
model is indicated by a white cross-point. Figure 4 also 
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Fig. 4. — UV to near-IR observed flux densities (erg s _1 cm -2 A -1 ) and best-fitting STARBURST99 model SEDs (using the Geneva 
(left) and Padova (right) evolutionary tracks) for the SSCs in SBS 0335-052. The blue data points (F220W, F330W, F435W, F550M) are 
included in the fitting process but the red points (F791W, F160W, F187N, F205W) are not. Grids of a (the standard deviation of the 
logarithmic residuals of the observed and model flux densities) are shown for models of various ages and extinctions. Black indicates better 
fitting models (lower values of a) and red indicates worse fitting models (higher values of a). A white cross indicates the best-fitting model. 
The dip in Ay at ~ 5 Myr is the result of a dramatic reddening of the intrinsic model colors (e.g. U — B) at this time. Lower values of 
Ay are preferred during this time since the intrinsic colors are already very red. 



shows that there is a clear excess in the /-band and at 
near-IR wavelengths compared to the model SEDs. We 
will return to this in §4.3. 

Estimates of the ages, masses, and extinctions of the 
SSCs are found by averaging these physical properties 
of the best-fitting model SEDs from the two STAR- 
BURST99 simulations and the results are listed in Table 
4. We note that the derived extinctions for SSCs 1 and 2, 
Ay ~ 0.5, are consistent with the extinction from the op- 
tical recombination line ratio found by Izotov et 



al. (1997). There is no significant difference between the 
goodness-of-fits using the two different STARBURST99 
simulations. The quoted errors reflect the variations in 
the results using the two simulations, although they rep- 
resent the minimum uncertainties since the models are 
somewhat degenerate in age and extinction (see Figure 
4). Despite this degeneracy, we still have confidence in 
the physical properties of the SSCs derived from SED 
fitting for the following reasons. First, the ages are con- 
sistent with those we derive from Ha equivalent widths 
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TABLE 4 
Properties from SED Fitting 



Source 


Age 


Mass 


A v 




(Myr) 


(1O 6 M ) 


(mag) 


SSC 1 


< 3 a 


1.0(0.1) 


0.5(0.1) 


SSC 2 


< 3 a 


1.1(0.1) 


0.5(0.1) 


SSC 3 


5.7(3.8) 


0.4(0.2) 


0.2(0.1) 


SSC 4 


< 3 a,b 


1.1(0.1) 


0.1(0.1) 


SSC 5 


1^8(1.1) 


1.8(0.2) 


0.0(0.1) 


SSC 6 


12.2(3.2) 


0.2(0.1) 


0.1(0.1) 



Note. — Estimates of the physical prop- 
erties of the SSCs in SBS 0335-052 derived 
from SED fitting. The values listed are 
averages from the best-fitting model SEDs 
using the two STARBURST99 simulations 
described in the text (there was no signif- 
icant difference in the fits using either of 
the model simulations). The quoted errors 
in parenthesis reflect the variations in the 
results using the two simulations. Also see 
Figure 4. 

a The SEDs are insensitive to ages < 3 
Myr. 

b This young age is probably incorrect 
since copious amounts of Ha and Pa emis- 
sion would be expected for such a young 
cluster and this is not observed. The age 
of SSC 4 inferred from the Ha equivalent 
width is ~ 12 Myr. 



in the following section (with the exception of SSC 4, see 
below). Also, as noted above, the Ay are in agreement 
with the work of Izotov et al. (1997), although we empha- 
size that caution must be applied when interpreting the 
meaning of extinction estimates of the embedded SSCs 
(1 and 2). We will address this issue in detail in §4.2.3. 

4.2. Nebular Emission from the SSCs 
4.2.1. Ages 

In addition to SED fitting to the stellar continuum, the 
nebular emission from a starburst region can also provide 
estimates of the physical properties of the host cluster. 
Age estimates can be obtained from the equivalent width 
of Ha emission since this quantity measures the ratio of 
ionizing flux (from massive, short-lived stars) to the total 
continuum flux density. This ratio is strongly dependent 
on age for stellar populations between approximately 3 
and 20 Myr old. The insensitivity before ~ 3 Myr arises 
because the most massive stars must begin to die in or- 
der for the ionizing flux (and Ha equivalent width) to 
decrease. Clusters older than ~ 20 Myr are not strong 
Ha emitters. 

In addition, if the dusty interstellar medium (ISM) sur- 
rounding the clusters is uniform and every Lyman con- 
tinuum photon ionizes a hydrogen atom, rather than get- 
ting absorbed by dust or escaping the cluster altogether, 
then Ha equivalent width will be independent of redden- 
ing (since the continuum is measured at the wavelength 
of Ha). If, however, the surrounding ISM is porous and a 
significant fraction of ionizing and non-ionizing photons 
escape the cluster, the equivalent width of Ha will be 
artificially low and only provide an upper limit on the 
age of the cluster. This is because we would detect all 
of the stellar continuum photons at 6563 A that escape, 



but the Ha flux would be reduced due to the leakage of 
ionizing photons (which do not actually ionize the gas). 

Age estimates from comparing the Ha equivalent 
widths of the SSCs with the STARBURST99 models are 
given in Table 5 and a significant age spread between 
the clusters is found (< 3 to ~ 15 Myr). These age 
estimates are consistent with those obtained from SED 
fitting (within the errors), with the exception of SSC 4. 
In this case, the age from SED fitting (< 3 Myr) is most 
likely incorrect since we would expect copious amounts of 
Ha and Paa emission from such a young cluster and this 
is not observed. The age of SSC 4 from the Ha equiv- 
alent width is ~ 12 Myr. It is not entirely clear why 
such a young age was found from SED fitting, although 
contamination by the diffuse emission surrounding SSC 
4 (see Figure 1) at short wavelengths may be responsible. 

4.2.2. Ionizing Luminosities 

The production rate of Lyman continuum photons, 
Qlvc? from the SSCs can be estimated from their Ha and 
Paa fluxes. From Equations 2 and 3 in Condon (1992), 
the ionizing photon rate as a function of H/3 luminosity, 
L H/ 3, is given by 

This equation provides only a lower limit since a frac- 
tion of the ionizing photons may be absorbed by dust, 
or escape the cluster altogether before ionizing hydro- 
gen atoms if the ISM is porous and clumpy. For Case 
B recombination, an electron temperature of 20,000 K 
(Izotov et al. 1997), and an electron density between 10 2 
and 10 4 cm -3 , the expected (reddening-free) flux ratios 
of Ha and Paa to R(3 are given by Condon (1992) and 
Osterbrock (1989): F Ha /F H0 = 2.73, F Paa /F up = 0.28. 
We use Equation 1 and the predicted line ratios with our 
observed (not extinction corrected) Ha and Paa fluxes 
to obtain estimates of the ionizing luminosities, Qh yc , 
produced by the SSCs in SBS 0335-052 (see Table 5). 

In addition to the Ha and Paa fluxes, we use new mea- 
surements of the (thermal) free- free radio flux densities at 
1.3 cm (Johnson, Hunt, & Reines in prep) to estimate the 
ionizing luminosities of SSCs 1 and 2 (the other clusters 
do not have significant detections at 1.3 cm). Johnson, 
Hunt, & Reines (in prep) find that the radio flux densi- 
ties at 1.3 cm are optically thin and equal to 0.17 mJy 
and 0.12 mJy for SSCs 1 and 2, respectively. We adopt 
a distance of 54 Mpc to calculate the radio luminosity 
densities, i^thermab an d follow Condon (1992) 

( -^V, thermal \ 

X VlO 27 erg s- 1 Hz- 1 J ' 

to estimate the ionizing luminosities of SSCs 1 and 2 (see 
Table 5). 

The ionizing luminosities of the SSCs should decrease 
as they age. Since we have established an age spread be- 
tween the different clusters, we can check this prediction. 
However, ionizing luminosity is also a function of mass 
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TABLE 5 

Properties from the Nebular Emission 





Age 
(Myr) 


nl .3cm 
^Lyc 

(10 49 s- 1 ) 


Qlyc 

(10 49 s- 1 ) 


^Ly*c 

(10 49 s- 1 ) 


fl .3cm 

T Q 


^Paa 

T Q 


T Q 


SSC 1* 


< 3.3 


> 3750 


> 2110 


> 1860 


0.76(0.16) 


0.43(0.09) 


0.38(0.08) 


SSC 2* 


< 3.4 


> 2650 


> 1350 


>1140 


0.51(0.11) 


0.26(0.06) 


0.22(0.05) 


SSC 3 


6^8(2.5) 




160(20) 


190(20) 




- 1 


~ 1 


SSC 4 


12.4(1.7) 




70(10) 


50(20) 




~ 1 


~ 1 


SSC 5 


15.1(2.3) 




50(20) 


40(30) 




- 1 


~ 1 


SSC 6 


13.9(1.9) 




20(10) 


10(10) 




- 1 


~ 1 



Note. — Age estimates of the SSCs from Ha equivalent widths, as well as measured ionizing 
fluxes (QLyc) an d recovered fractions of ionizing fluxes (/q) from Ha, Paa and 1.3 cm radio flux 
densities (Johnson, Hunt, & Reines in prep). A three-dot ellipsis indicates sources that do not 
have detections at 1.3 cm. 

*The age estimates for these young clusters are upper limits since Ha equivalent width is insen- 
sitive to ages < 3 Myr and also possibly due to the leakage of ionizing photons through a porous 
ISM (see discussion in the text). The ionizing fluxes for SSCs 1 and 2 are lower limits because they 
do not reflect ionizing photons that either escape through a porous ISM or are absorbed by dust. 
Correcting the Q\, yc by the /q's imply SSCs 1 and 2 each have ionizing fluxes of ~ 5 x 10 52 s -1 . 




5 10 15 20 

Age [Myr] 



Fig. 5. — Scaled ionizing fluxes versus age (from Ha equivalent 
width). For SSCs 1 and 2, the ionized gas emission underpredicts 
the expected values from the non-ionizing stellar continuum. This 
suggests the ISM is porous and stellar continuum photons (both 
ionizing and non-ionizing) are leaking out of the immediate vicinity 
of the clusters. See the text for a full discussion. 



(in addition to age) so we need to scale the measured 
QLyc values to a common mass and remove the mass de- 
pendence. To achieve this, we divide the measured Q\, yc 
values by the masses derived from SED fitting (Table 4) 
and multiply by 10 6 M to scale the ionizing luminosi- 
ties to those expected for a 10 6 M cluster. Figure 5 
shows the scaled ionizing luminosities versus cluster age 
(derived from Ha equivalent widths). The two STAR- 
BURST99 models are also shown for comparison. 

From Figure 5, we see that the ionizing luminosities of 
SSCs 3-6, derived from their Ha and Paa fluxes, are con- 
sistent with the model predictions. However, the ioniz- 
ing luminosities of SSCs 1 and 2 (ages < 3 Myr), derived 
from their Ha and Paa fluxes, are well below the ex- 
pected values. In other words, these recombination lines 
are only sampling a fraction of the ionizing photons rela- 
tive to the number expected from the non-ionizing stellar 
continuum inferred from the UV and optical photometry 



(recall that we have normalized the Qhy C values by the 
masses derived from SED fitting to the UV and optical 
photometry) . Equivalently, the masses inferred from the 
recombination line fluxes are less than those derived from 
SED fitting. 

Interestingly, the optically-thin, free- free radio emis- 
sion from SSCs 1 and 2 (Johnson, Hunt, & Reines in 
prep) also underpredicts the ionizing luminosities (see 
Figure 5). This may suggest the ISM surrounding the 
young SSCs is clumpy and porous, with virtually unob- 
scured lines-of-sight to the cluster stars. Clearly, some 
fraction of stellar continuum photons escape the clus- 
ters through low-extinction regions: we observe copious 
amounts of UV and optical photons and the Ay derived 
from SED fitting is a modest ~ 0.5 mag for SSCs 1 and 
2. The lower than expected QLyc suggest ionizing pho- 
tons also escape the immediate vicinity of the clusters 
(without ionizing hydrogen), resulting in the loss of Ha, 
Paa and radio emission. This in turn, implies the Ha 
equivalent widths are lower limits and the ages of SSCs 
1 and 2 are < 3 Myr. 

We can calculate the recovered fractions of the ionizing 
luminosities from SSCs 1 and 2 by comparing the data 
points to the models in Figure 5. The fraction of ionizing 
luminosity recovered by Ha, /g a , is ~ 0.38 for SSC 1 
and r^j 0.22 for SSC 2. These values are consistent with 
previous estimates that, in the optical, only about 1/4 of 
the total ionized gas emission is observed for the region 
containing SSCs 1 and 2 (Hunt et al. 2001). From Paa, 
yPaa = q 43 ( SSC ^ and o 26 (SSC 2). Using the radio- 
derived ionizing luminosities, we find /g 3cm equals ~ 
0.76 and ~ 0.51 for SSCs 1 and 2, respectively. The 
increase of /q with wavelength naturally follows from 
differential absorption (i.e. reddening) of the ionized gas 
as probed by the Ha, Paa and radio emission. 

The measured ionizing luminosities and the recovered 
fractions of the intrinsic ionizing luminosities are sum- 
marized in Table 5. Correcting the measured ionizing 
luminosities (i.e. dividing the Qlvc's by /q' s ) implies 
the intrinsic ionizing luminosities of SSCs 1 and 2 are 
each - 5 x 10 52 s _1 , or the equivalent of - 5000 07.5 V 
stars (Vacca et al. 1996). The inferred masses of these 
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SSCs are - 1O 6 M (see Table 4). 

4.2.3. Extinctions of the Embedded SSCs 

The youngest SSCs in SBS 0335-052 (1 and 2) are still 
embedded in their natal birth material: they are the 
strongest emitters of Ha and Paa and the only clusters 
in the galaxy detected in the radio via free-free emission 
(Johnson, Hunt, & Reines in prep). Various studies have 
provided extinction estimates from recombination lines 
of the region encompassing these two clusters (previous 
studies were not able to resolve the two sources) and the 
derived values appear to be discrepant: Ay =0.55 from 
Ha/H/3 (Izotov et al. 1997), Ay = 0.73 from H^/Br 7 
(Vanzi et al. 2000), A v = 1.45 from ftp /Bra (Hunt et 
al. 2001), and A v = 12.1 from Br7/Bra (Hunt et al. 
2001). Mid-IR observations complicate the situation fur- 
ther. Thuan et al. (1999) found that the mid-IR SED is 
well-described by a modified blackbody spectrum extin- 
guished by a screen of dust with Ay ~ 20 mag. A qual- 
itatively similar result was found by Plante & Sauvage 
(2002) and Hunt et al. (2005). These authors model the 
observed mid-IR SED using DUSTY (Ivezic & Elitzur 
1997; Ivezic et al. 1999), a program that solves the radi- 
ation transfer equations in a spherical environment, and 
conclude that the embedded starburst suffers ~ 30 mag- 
nitudes of visual extinction. Houck et al. (2004) also find 
a large extinction, Ay ~ 15, from the 9.7 /im silicate 
absorption feature assuming a screen model and a black- 
body background source. In this section, we attempt to 
reconcile the apparently discrepant extinction estimates 
assuming all of the emission (at all wavelengths) is asso- 
ciated with the optically visible SSCs 1 and 2. 

We begin by estimating the extinctions of the recom- 
bination lines R(3 (Izotov et al. 1997), Ha (this work), 
Paa (Thompson et al. 2006, and this work), Br7 (Vanzi 
et al. 2000), and Bra (Hunt et al. 2001) relative to the 
(thermal) free-free radio emission from this region (John- 
son, Hunt, & Reines in prep). The radio emission does 
not suffer from extinction per se, although again we note 
that the far-UV photons can be absorbed by dust or es- 
cape the boundaries of the clusters before ionizing the 
gas. In either of these scenarios, the ionizing flux probed 
by the optical and IR lines would also be reduced by the 
same amount and the relative extinction between these 
lines and the radio would not be affected. To obtain the 
extinction estimates, we compare the measured fluxes to 
the expected values predicted by the free-free radio emis- 
sion at 1.3 cm. Following Condon (1992), the predicted 
Ha flux is given by 

(5?^)~ 08x1 »- ,2 (#k)" 58 < 3 » 

/ V \ 0,1 ( ^thermal A 

x vghzJ V m Jy / ' 

where ^thermal is the thermal free-free radio flux den- 
sity (= 0.56 mJy). To obtain the predicted fluxes for 
the other recombination lines, we use the theoretical 
(reddening-free) line ratios appropriate for SBS 0335-052 
(Case B, T=20000 K) given by Condon (1992) and Os- 
terbrock (1989). The observed recombination line fluxes 
and the radio-predicted values are listed in Table 6. We 
give the flux ratios (observed to predicted) and the corre- 



TABLE 6 

Recombination Lines from the Region SSC 1+SSC 2 



Line Rest A Observed Predicted Flux A\ 





w 


Flux a 


Flux a 


Ratio 


(mag) 




0.486 


6.1 


14.3 


0.4 


0.9 


Ra 


0.656 


16.6 


40.8 


0.4 


1.0 


Paa 


1.875 


1.77 


4.05 


0.4 


0.9 


Br 7 


2.166 


0.15 


0.34 


0.4 


0.9 


Bra 


4.052 


0.90 


0.88 


1.0 


0.0 



Note. — Observed recombination line fluxes for H/3 
(Izotov et al. 1997), Ha (this work), Paa (Thompson et 
al. 2006, and this work), Brj (Vanzi et al. 2000), and 
Bra (Hunt et al. 2001) were measured in a ~ l'/5 x l'.'O 
aperture. The predicted values are from the free-free radio 
flux density at 1.3 cm (Johnson, Hunt, & Reines in prep). 

a Units are 10 -14 erg s _1 cm -2 . 



sponding extinctions of these lines at their native wave- 
lengths, A\. 

It is quite intriguing that, relative to the free- free radio 
emission, the calculated extinctions of the Ha, Paa, 
and Br7 emission are all ~ 1 mag (at their native wave- 
lengths between ~ 0.5 — 2 /im). The similar estimates 
of A\ (or equivalently, the observed to predicted flux ra- 
tios) for these lines appear to indicate that there is not a 
significant amount of reddening (i.e. differential absorp- 
tion with short wavelengths suffering more absorption) 
of the ionized gas between ~ 0.5 — 2 /im. 

The value A\ ~ 1 then gives the total amount of ab- 
sorption of these lines compared to the radio (i.e. the 
ratio of the ionized gas probed by the optical/near-IR 
lines to the total amount of ionized gas inferred from 
the free- free radio emission). If the ISM is clumpy as we 
have previously suggested (§4.2.2) and no optical or near- 
IR light gets through dense dust clumps (i.e. the light 
is completely absorbed), the ratio of the optical (and 
near-IR) emission to the radio emission is equivalent to 
the ratio of the amount of flux that emerges through 
the "holes" (optical/near-IR) to the amount of flux that 
emerges through both the holes and the clumps (radio). 
This ratio equals 0.4 for Ha, Paa, and Br7« In other 
words, ~ 40% of the dust cocoon is void of clumps (low 
extinction regions) and ~ 60% is covered by clumps of 
dust that totally absorb light shortward of ~ 2 /im (high 
extinction regions). In effect, there is a clumpy dust cov- 
ering factor of ~ 60%. 

These values are strikingly similar to our estimates of 
the fraction of ionizing luminosity recovered by the radio 
emission (~ 64% on average for SSCs 1 and 2) and the 
fraction of ionizing luminosity that escapes the clusters 
undetected (~ 37%, see §4.2.2). This may suggest that 
the dust clumps and the neutral gas clumps are coinci- 
dent, although this is not necessarily the case. 

The picture outlined above is illustrated schematically 
in Figure 6. The large extinctions (Ay > 15 — 30) derived 
from mid-IR observations (Thuan et al. 1999; Plante & 
Sauvage 2002; Houck et al. 2004; Hunt et al. 2005) would 
be a natural consequence if the mid-IR emission orig- 
inates from the dense dust clumps and not a uniform 
screen. In our proposed scenario, the dust would be 
heated by the strong UV stellar continuum impinging on 
the clumps. This would help explain the large extinction 
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Fig. 6. — A cartoon illustrating the effects of a clumpy dust 
structure surrounding the young SSCs (1 and 2) in SBS 0335-052. 
Short- wavelength photons (blue) with A < 2 /im (i.e. stellar con- 
tinuum, Hf3, Ha, Paai, and Br7) are completely absorbed by dense 
dust clumps harbouring large grains and only emerge through the 
diffuse interclump regions, accounting for the low measured extinc- 
tions in the optical/near-IR. Longer- wavelength photons (red) with 
A > 4 /im (i.e. Bra: and 1.3 cm) can penetrate the dense clumps 
as well as the diffuse regions. See §4.2.3 for a full discussion. 

derived from the mid-IR spectral energy distributions. 
Also, if the extinctions inferred from the mid-IR describe 
absorption from dust clumps rather than a screen of dust, 
they may actually underestimate the absorption because 
of beam dilution effects. In addition, diffuse dust in the 
interclump regions could provide the modest reddening 
measured at optical wavelengths (since the optical light 
only emerges through the low extinction regions). The 
scenario described here is similar to the one proposed by 
Gordon et al. (1997) in which a dust-free star cluster is 
surrounded by a clumpy shell of dust and gas. 

Moreover, the extinction curve in SBS 0335-052 ap- 
pears to be anomalous. It is roughly constant up to 2 /im 
with A\ ~ 1, then drops to A\ ~ at 4 /im. This be- 
havior is qualitatively similar to extinction curves along 
lines of sight with high gas densities. Such regions have 
typically natter extinction curves, thought to be due to 
larger dust grains than in the diffuse ISM (e.g., Kim et 
al. 1994). In the case of SBS 0335-052, the grains may be 
as large as ~ 1 — 2 /im; however the silicate absorption 
feature at 9.7 /im (e.g., Houck et al. 2004) makes grains 
larger than this rather unlikely, because in that case the 
absorption feature would be absent (Maiolino & Natta 
2002). 

4.3. I -band and Near-IR Excesses 

The SSCs in SBS 0335-052 exhibit excesses in the I- 
band and at near-IR wavelengths with respect to model 
SEDs (Figure 4). Excesses in the F791W (~ J), F160W 
(~ H), F187N (~ 1.9 /im continuum), and F205W (~ K) 
filters are given in Table 7, where an excess in mag- 
nitudes is defined as 2.5xlog(/ observed // mode i) (/observed 
and /model are the observed and best-fitting model flux 
densities (§4.1), respectively). In this section, we inves- 
tigate the nature of these excesses. We emphasize that 
the /-band (~ 0.8 /im) and near-IR (~ 1.6 — 2.1 /im) 
excesses likely have different origins. 



TABLE 7 
Near-IR Excesses of the SSCs 



Source 


F791W 


F160W 


F187N 


F205W 


SSC 1 


0.7 


1.2 


2.0 


2.7 


SSC 2 


0.5 


1.3 


1.8 


2.6 


SSC 3 


0.6 


1.5 


1.6 


2.0 


SSC 4 


0.5 


0.3 


0.8 


1.1 


SSC 5 


0.5 


0.4 


0.9 


1.1 


SSC 6 


0.6 


0.6 


0.8 


1.3 



Note. — Near-IR excesses (in magnitudes) of 
the SSCs with respect to the best-fitting model 
SEDs (to the UV and optical data). An excess 
is defined as 2.5xlog(/ bserved//model), where 
/observed an d /model are the observed and model 
flux densities, respectively. The values listed are 
averages from the best-fitting model SEDs using 
the two STARBURST99 simulations described in 
the text (there was no significant difference in the 
fits using either of the model simulations). 



4.3.1. Excess in the I -band: Extended Red Emission? 

The /-band excess observed in the clusters of SBS 
0335-052 ranges from ~0.5-0.7 mag, similar to the excess 
found in the young massive clusters of NGC 4449 8 (~ 0.5 
mag on average; Reines et al. 2008). The primary origin 
of the /-band excess in these two systems is quite possibly 
the same. An extensive discussion of possible origins of 
the /-band excess observed in the young clusters of NGC 
4449 is given in Reines et al. (2008) and the favored hy- 
pothesis is Extended Red Emission (ERE) (for a review 
of ERE, see Witt & Vijh 2004). ERE is a photolumi- 
nescent process in which an interstellar particle absorbs 
a UV photon and then ultimately emits a photon in the 
~ 6000 — 9000 A wavelength range. ERE is typically ob- 
served as a broad (~ 1000 A) emission feature and it has 
been detected in many astrophysical environments where 
both dust and UV photons are present, including the 
Orion Nebula (Perrin & Sivan 1992) and the 30 Doradus 
Nebula (Darbon et al. 1998). ERE provides a favorable 
hypothesis for the /-band (~ 0.8 /im) excess observed 
in the SSCs of SBS 0335-052, but it cannot explain the 
excesses seen at near-IR wavelengths (~ 1.6 — 2.1 /im). 

4.3.2. Excess in the Near-IR: Hot Dust in the Younger 
SSCs and Red Supergiants in the Older SSCs 

A near-IR/optical color-color diagram can provide in- 
sight into the origin of the near-IR (~ 1.6 — 2.1 /im) 
excesses observed in the SSCs of SBS 0335-052. Fig- 
ure 7 shows a plot of [F550M]-[F205W] (~ V — K) 
versus [F550M]-[F160W] (~ V - H) with the SSCs la- 
beled. The line [F160W]-[F205W]=0 is also shown, ap- 
proximating the color of main sequence stars through- 
out their lives and consistent with the STARBURST99 
evolutionary models. In addition, we plot the line 
[F160W] — [F205W]=0.5, a characteristic red supergiant 
(RSG) color (McCrady et al. 2003; Pickles 1998) 9 . Note 

8 We note that the /-band observations of SBS 0335-052 and 
NGC 4449 come from different instruments and slightly different 
filters (WFPC2/F791W for this work and ACS/F814W for NGC 
4449). 

9 McCrady et al. (2003) calculate a synthetic [F160W]-[F222M] 
color of 0.52 for a M2 I star in the Pickles (1998) library. 
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Fig. 7. [F550M]-[F205W] (~ V-K) versus [F550M]-[F160W] 
(~ V — H) color-color diagram. The F550M magnitudes are dom- 
inated by the young hot stars in the clusters, while the F160W 
and F205W magnitudes of all of the SSCs are much brighter than 
predicted from the STARBURST99 models (Table 7). SSCs 3-6 
(ages > 7 Myr) have [F160W]-[F205W] colors consistent with red 
supergiants, while SSCs 1 and 2 (ages < 3 Myr) have much redder 
colors indicating another origin (see Figure 8). 

that in this color-color diagram, the F550M magnitudes 
are dominated by the hot, unevolved cluster stars. 

Figure 7 illustrates that SSCs 3-6 have 
[F160W]-[F205W] colors consistent with RSGs. 
This is in agreement with other studies showing these 
evolved stars are the dominant source of near-IR light 
in SSCs older than ~ 7 Myr (e.g. McCrady et al. 
2003). We should also point out that RSGs could also 
contribute to the observed /-band excess in these older 
clusters. The two youngest SSCs in SBS 0335-052 (1 
and 2) have extremely red [F160W] — [F205W] colors 
that are quantitatively different from the other clusters: 
~ 1.5 and ~ 1.3, respectively (Table 2). The origin of 
these red colors is not consistent with RSGs and another 
mechanism must be responsible. This is not surprising 
since RSGs should not have had enough time to evolve 
in these < 3 Myr old clusters. 

Emission lines have been observed in the H and K 
spectra of SBS 0335-052 (Vanzi et al. 2000) in the region 
containing SSCs 1 and 2, however the line contribution 
to the broad-band photometry here is negligible. In ad- 
dition, the near-IR excesses observed in SSCs 1 and 2 
(Figure 4) are so large (~ 1.3 mag in F160W, ~ 1.9 
mag in F187N, and - 2.7 mag in F205W; Table 7), and 
an excess is present in the narrow-band continuum filter 
F187N, that continuum emission must be the dominant 
cause. 

Although ionized gas continuum emission may be im- 
portant (see e.g. Vanzi et al. 2000; Hunt et al. 2001), 
we propose that hot dust is the main origin of the red 
[F160W]-[F205W] colors for SSCs 1 and 2. To inves- 
tigate this hypothesis further, we add a modified black- 
body with a dust opacity curve given by Whitney et al. 
(2003) and Kim et al. (1994) to the best-fitting model 
SEDs (to the UV and optical data) presented in §4.1. 
The results are positive and are shown in Figure 8. Tem- 
peratures of 800 and 850 K for the hot dust components 
in SSCs 1 and 2 respectively, are effective at reproduc- 
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Fig. 8. — UV to near-IR spectral energy distributions for SSCs 
1 and 2. For each cluster, we add a hot dust component (modified 
blackbody) to the best-fitting STARBURST99 model (to the UV 
and optical data) found in §4.1. Temperatures equal to 800 K 
(SSC 1) and 850 K (SSC 2) are effective at reproducing the near- 
IR excesses. 

ing the observed near-IR SEDs. Vanzi et al. (2000) also 
concluded that hot dust was necessary to explain the red 
H — K colors they observed in SBS 0335-052. 

5. EVIDENCE FOR SUCCESSIVE CLUSTER FORMATION 
IN SBS 0335-052 

We have already established a significant age spread 
between the SSCs in SBS 0335-052 (§4.2.1). In this sec- 
tion, we investigate the possibility of successive cluster 
formation in the galaxy. Thuan et al. (1997) demon- 
strated a systematic increase in the [F569W] — [F791W] 
color away from SSC 1. They attributed this trend to a 
combination of differential extinction by dust and evolu- 
tionary effects arising from sequential (self) propagating 
star formation. Thompson et al. (2006) also suggest that 
each star-formation event has positive feedback on sub- 
sequent star formation. 

In the top of Figure 9, we plot [F550M]-[F791W] 
color versus projected distance from SSC 1, similar to 
Figure 2 in Thuan et al. (1997). We use the same /- 
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shocks), although this is a possibility. Alternatively, the 
SSCs in SBS 0335-052 may have been sequentially trig- 
gered by a large-scale disturbance sweeping through the 
galaxy from roughly north to south, the origin of which 
is presently unknown, but could be related to an interac- 
tion with neighboring galaxies SBS 0335-052W and NGC 
1376 (Pustilnik et al. 2001). 

6. CONCLUSIONS 

We have presented a multi- wavelength study of the su- 
per star clusters in SBS 0335-052, the lowest metallicity 
galaxy with a star formation rate > 1 M yr _1 known. 
We use new near-IR and archival optical HST observa- 
tions as well as free-free radio continuum measurements 
(Johnson, Hunt, & Reines in prep) of the galaxy to probe 
the stellar populations and the gas and dust components 
of the clusters. The main results of our work are sum- 
marized below: 
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Fig. 9. Top: [F550M]-[F791W] (V - I) color versus projected 
distance from SSC 1. Observed colors and colors corrected for the 
/-band excess (Table 7) are shown. We do not observe a correlation 
between V — I and distance, as Thuan et al. (1997) did. This is 
due to the different V-band filters used: the F569W filter used by 
Thuan et al. (1997) is highly contaminated by gaseous emission 
lines, whereas the F550M filter used here is not. Bottom: Age 
(from Ha equivalent width) versus projected distance from SSC 1. 
There is a strong correlation indicating successive cluster formation 
in SBS 0335-052 from the northern (SSC 6) to southern (SSC 1) 
end of the galaxy. 



band (F791W) filter here, but a different (line-free) V- 
band filter (F550M). The F569W filter used by Thuan 
et al. (1997) contains significant contamination from 
gaseous emission lines, whereas the F550M filter does 
not. We do not observe the same strong trend of redder 
[F550M]-[F791W] colors with distance from SSC 1. 

We do, however, observe a strong correlation between 
cluster age (from Ha equivalent width) and distance from 
SSC 1 in the bottom of Figure 9. This is clear evidence 
for successive cluster formation from the northern end 
of SBS 0335-052 (SSC 6) to the southern end of the 
galaxy (SSC 1): we calculate a "propagation" speed of 
— 35 pc Myr -1 w 35 km s _1 . It should be emphasized 
that successive cluster formation does not necessarily im- 
ply self-propagating cluster formation (i.e. SSCs trig- 
gering the formation of more SSCs through supernovae 



1. There is a significant age spread between the SSCs 
in SBS 0335-052 that is correlated with position in 
the galaxy. The ages are in the range < 3 Myr 
to ~ 15 Myr with the youngest clusters located 
in the south and the oldest clusters in the north. 
A large-scale disturbance, with a velocity of — 35 
km s _1 , appears to have triggered the successive 
cluster formation that we observe. 

2. We find evidence for a porous and clumpy ISM sur- 
rounding the youngest SSCs in SBS 0335-052 (1 
and 2). The measured ionizing luminosities from 
Ha, Paa, and optically-thin, free-free radio emis- 
sion (Johnson, Hunt, & Reines in prep) are lower 
than expected compared to the optical SEDs, sug- 
gesting a fraction of ionizing photons from the stel- 
lar continuum are escaping the strict confines of 
the clusters before ionizing hydrogen, that would 
have otherwise contributed to the measured ionized 
gas emission (i.e. Ha, Paa, and radio). The Ay's 
of SSCs 1 and 2 derived from SED fitting (~ 0.5 
mag) also indicate the existence of low-extinction 
regions that provide relatively unobscured lines-of- 
sight into the embedded clusters. The corrected, 
intrinsic ionizing luminosities of SSCs 1 and 2 are 
each - 5 x 10 52 s" 1 (the equivalent of - 5000 07.5 
V stars) and they each have total stellar masses of 
- 1O 6 M . 

3. An /-band (F791W, ~ 0.8 /im) excess with re- 
spect to model SEDs is observed for all of the 
SSCs in SBS 0335-052. The /-band magnitudes 
are ~ 0.5 — 0.7 mag brighter than the model predic- 
tions. A similar result was found in the young mas- 
sive clusters in NGC 4449 (Reines et al. 2008) and 
attributed to a photoluminescent process known as 
Extended Red Emission. We hypothesize that the 
same mechanism dominates the /-band excess in 
SBS 0335-052. 

4. All of the SSCs have red near-IR (-1.6-2.1 /im) 
colors and large near-IR excesses with respect to 
model SEDs fit to the optical photometry. The 
red near-IR colors, however, are quantitatively dif- 
ferent (redder) in the youngest SSCs (1 and 2) 



SSCs in SBS 0335-052 



13 



compared to the older SSCs (3-6), clearly indicat- 
ing a different origin for the red colors. We have 
demonstrated that the near-IR light from the older 
clusters is dominated by evolved red supergiants, 
whereas the near-IR light from the youngest clus- 
ters is dominated by hot dust emission. The young 
SSCs 1 and 2 are also the only radio-detected clus- 
ters in the galaxy (Johnson, Hunt, & Reines in 
prep). 

5. We have proposed a scenario that can account 
for the apparently discrepant extinction estimates 
found in the literature for the starburst region en- 
compassing the youngest SSCs (1 and 2) in SBS 
0335-052. Our picture is consistent with all of the 
emission (at all wavelengths) being associated with 
the optically visible SSCs 1 and 2, and is based on 
our calculated absorptions of the optical and IR re- 
combination lines relative to the (thermal) free-free 
radio emission (Johnson, Hunt, & Reines in prep). 
We find that an ISM containing dense dust clumps 
harbouring large grains (~ 1 — 2 /im), with a dust 
clump covering factor of ~ 60%, can account for 
the absorption of the optical/near-IR recombina- 
tion lines relative to the radio, as well as the large 
extinctions derived from previous mid-IR observa- 
tions. Interclump regions containing diffuse dust 
can account for the low measured extinctions de- 
rived from optical/near-IR observations. 

The main goal of this case study was to investigate 
the formation of massive star clusters in an extremely 
low metallicity environment that might be analogous to 
conditions in which the now ancient globular clusters 
were born. In general, the formation of massive star 
clusters in SBS 0335-052 appears to be similar to that 
found in higher metallicity counterparts, e.g. Henize 2- 
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10 and Haro 3 (Johnson & Kobulnicky 2003; Johnson 
et al. 2004). These galaxies are all blue compact dwarfs 
(BCDs), and each of them currently hosts radio-detected 
massive star clusters with ages < 3 Myr and masses 
of rsj 1O 6 M , all of which are consistent with evolving 
into globular clusters in several billion years. Perhaps 
even the low metal abundance in SBS 0335-052 is insuf- 
ficiently extreme to reflect primordial star formation; in 
fact, models predict that the transition from metal- free 
massive Population III stars to "normal" (less massive) 
Population II occurs at Z/Zq < 10 -5 (Mackey et al. 
2003; Bromm & Larson 2004). On the other hand, large 
metallicity fluctuations because of radiative feedback ef- 
fects may enable "hidden" Population III star forma- 
tion at slightly higher metal abundances (Tornatore et 
al. 2007). Clearly more detailed studies need to be car- 
ried out on metal-poor star-forming galaxies in the Local 
Universe, and comparative samples need to be obtained 
in order to investigate interactions and trends among star 
cluster and dust formation, environment, and metal en- 
richment. 
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